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Abstract 

The crack growth properties of several sealing glasses were measured using 
constant stress rate testing in 2% and 95% RH (relative humidity). Crack growth 
parameters measured in high humidity are systematically smaller (n and B) than those 
measured in low humidity, and velocities for dry environments are ~ 1 OOx lower than for 
wet environments. The crack velocity is very sensitivity to small changes in RH at low 
RH. Confidence intervals on parameters that were estimated from propagation of errors 
were comparable to those from Monte Carlo simulation. 

Introduction 

Sealing glasses are used in components such as electrical feed through connectors. 
The glass seals and electrically insulates the connector, and thus fracture of the brittle 
seal is a concern. In applications such as the Space Shuttle Environmental Cut Off 
(ECO) system, the connector seals are subjected to differential pressures at cryogenic 
temperatures and seal failure can create leakage of dangerous liquids and/or gasses. 

The crack growth properties of several sealing glasses were measured to compare 
glasses and to help perform life prediction and reliability of components such as feed 
through connectors. Strength based measurements, which are convenient, were used to 
generate the data. However, because the statistical scatter in parameters derived from 
strength data can be very large, the statistical significance of the estimates was checked 
by estimation of confidence intervals on the parameters via propagation of errors (POE) 
and Monte Carlo methods. The large scatter is a result of strength not being a material 
property for glasses, but a function of the fracture toughness and worst flaw size present 
from a variety of sources. Ideally, parameter estimation and design of brittle materials 
should be done on a fracture mechanics basis (e.g. NASGRO [1]) rather than a strength 
basis because strength is a function of the highly variable flaw size and relatively 
consistent fracture toughness. 

Although fracture mechanics specimens with large cracks, such as the double- 
torsion specimen, can be used to measure crack growth with less scatter, the results are 
complicated by R-curve effects in coarse grain materials such as ZnSe [2] and diffusion 
rate effects when the crack size is large relative to that in real components. Strength 
based testing can be made more akin to fracture mechanics methods by placing a small 
precrack, such as an indentation, in specimens and thereby reduce scatter, yet test cracks 



on the order of those encountered in applications. Comparison of sealing glass 
parameters from strength and fracture mechanics methods is left to future study. 

In order to cover the range of environments to which components with sealing 
glasses are exposed, RH (relative humidity) of ~2% and 95% were considered. To 
expedite the work, constant stress rate testing of flexure specimens was used. The data 
was analyzed by linear regression of (1) the individual data points, (2) the median values, 
and (3) the average values. 


Materials 

The materials tested 1 were Corning 0120 sealing glass, Electro-Glass 2164, 
Schott 8330 borofloat glass, and Schott S-8070 SB glass-ceramic. In addition, the 
fracture toughness of several other glasses was measured for comparison: soda-lime 
silicate, S8061 sealing glass, and an antimony (Sb) doped glass. With the exception of 
the as-molded 2164 glass, the test specimens were prepared by diamond grinding in 
conformance with ASTM Cl 161 [3]. For the 2164 glass specimens for crack growth 
testing, the tensile surface was preserved in the as-molded condition. 

Experimental Procedure 

The elastic modulus of 2164 glass was determined at 20°C by impulse excitation 
of vibration in accordance with ASTM C 1259 [4], The mean and standard deviation of 
seven specimens in the as-molded and ground conditions were 62.0 ± 1.2 and 63.8 ± 0.5 
GPa respectively. 

Fracture strength as a function of stress rate was measured at 20°C by using four 
point flexure of ASTM Cl 161 size B specimens [3] at rates ranging from 10' to 1 0 
MPa/s in relative humidity ranging from ~2% to 95%. Typically, six stress rates were 
applied with at least five specimens per rate. For the purposes of parameter analysis, the 
inert strength was detennined by testing at low RH (<2%) with a stress rate greater than 
or equal to 1000 MPa/s. This results in failure in a fraction of a second. 

Fracture toughness was measured by using chevron-notch flexure specimens [5] 
in laboratory ambient (-30% RH) air and dry nitrogen. Test specimen stability was 
monitored via a strain gage placed on the compressive face of the specimen [6]. 

Data Analysis 


The power law formulation: 

v - d i- AK ’- A ' l t r <•> 

was applied in the data analysis, where v, a, and t are crack velocity, crack size, and time, 
respectively. A and n are the material/cnvironment dependent SCG (slow crack growth) 
parameters, and Ki and K IC are, respectively, the Mode / stress intensity factor and the 
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critical stress intensity factor or fracture toughness of the material under Mode / loading. 
For constant stress rate testing based on the power law formulation, the fracture 
strength, cy, is expressed as a function of stress rate as [7] 

a f = [B{n + lyrT 2 ^" 1 (2) 


where & is the applied stress rate, cr, is the inert strength, and B is a parameter associated 
with A, n, fracture toughness, crack geometry and loading configuration. The SCG 
parameter n can be determined from a plot of log Of as a function of log & with Eq. 2 
written as 


where 


log <J f = log & + log D 

n + 1 
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Once the slope a and intercept ft are estimated by linear regression of Eq. 3, the 
parameters n, D, B and A, and their standard deviation SD n , etc., were estimated from [8] 

n = — — 1 (5) 
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where Q = a- J3lnl0 + In cr i 

and Cov{a,f$)= -SD 2 a {logd) (15) 

where log d is the mean of the log of the applied stressing rates, Y is the geometry 
correction factor for the stress intensity factor, and the standard deviation associated with 
the inert strength (SDi nai ) is calculated in logarithmic space. Probability limits on the 
parameters B and A can be calculated from: 

B Upr „. = EXF{ln B ± t(SD,„ B )] and A Uppp , = EXf[l„ A ± t(SD „ J] (16) 

Lower Lower 

by using Student’s t distribution for the DOF and probability level desired. If the DOF 
(degrees-of-freedom) is greater than ~40, then 

B Um , =EXP[lnB±e(SD„ B )\ and A Upr „. =EXr\lnA±l(SD„ A )} (17) 

Lower Lower 



where t is the number of standard deviations corresponding to the probability level 
desired. The DOF, $ is given by 





where tp a is the DOF in inert strength (number of inert strength tests minus one) and (f) u p 

is the DOF in regression (number of constant stress rate tests minus two). 

Three approaches were used to estimate the slope and intercept of eq. 3: linear 
regression of (1) the individual data points; (2) the median values; and (3) the average 
values. In addition to the approaches described, the fits were performed over several 
stress rate ranges to determine the sensitivity to data range. 

Results 

Fracture Toughness 

Examples of load-backface strain curves for laboratory air and dry N 2 are shown 
in Figure 1 for the Electro-Glass 2164. Stable fracture was exhibited in both 
environments; however, less stability was exhibited in dry N?. Fracture toughness in dry 
N 2 ranged from 0.72 to 0.80 MPaVm for the glasses, as summarized in Table 1. Testing 
in air reduced the measured fracture toughness by -10%. For the purposes of crack 
growth parameter estimation, the fracture toughness the 0120, 8070 and 8330 materials, 
which were not measured, were assumed to be 0.75 ± 0.04 MPaVm. 
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Figure 1 Load as a function of backface strain for Electro-glass 2164 chevron-notched flexure specimens 

in diy nitrogen and laboratory’ air. 



Table 1 Fracture toughness (MPaVm) of several glasses 


Material 

Environment 


Air (%RH/°F ) 

Dry N 2 

S8061 

0.64 ± 0.01 (23/73) 

0.72 ± 0.02 

2164 

0.61 ±0.05 (32/73) 

0.74 ±0.03 

Soda lime silicate 

0.75 ±0.04 (35/73) 

0.80 ±0.01 

Sb-doped 

0.72 ± 0.002 (23/73) 

0.76 ±0.01 


Inert and Time-dependant Strength 

The strength as a function of stress rate is plotted in Figures 2-5. The large degree 
of scatter, particularly at low RH, is indicative of the difficulty in characterizing and 
designing glasses and dense optical materials with the inherent flaw population and 
strength measurements: random and spurious damage make the distribution ever 
changing and difficult to characterize, regardless of Weibull statistics. In this testing, the 
effect of scatter on slow crack growth was mitigated partially by the large range of stress 
rates used (>4 orders of magnitude). All the materials, except the 8070 SB glass-ceramic, 
exhibit a strength increase from ~50 MPa in 95% RH to —150 MPa in 2% RH as the 
stress rate is increased from 0.001 MPa/s to 1000 MPa/s, implying a similar combination 
of flaw size distribution and fracture toughness. As the fracture toughness values are 
likely similar (Table I), the implication is a similar flaw size distribution. 

The slow crack growth parameters as estimated from equations (5) to (17) are 
summarized in Table 2. 
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Figure 2 Strength of 0120 glass in 2% and 95% relative humidity. 
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Figure 3 Strength of 2164 in 2% and 95% relative humidity. 
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Figure 4 Strength of 8070 SB glass-ceramic in <2%> and 95% relative humidity. 
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Figure 5 Strength of 8330 borofloat in 1% and 95% relative humidity. 


Table 2 Summary of slow crack growth (SCG) parameters for glasses 


Regression of 
Individual Points 

11 

B 

MPa 2 s 

B-95%, 

A 

m/s 

(MPaVm)" n 

A+9S% r 

m/s 

(MPaVm)' n 

# of 
SCG 
Tests 

0120, 95% 

17.0 ±3.1 

0.6 

1.8 X 10" 4 

4.4 x 10 +0 

3.9 x 10 +4 

36 

0120, 2% 

23.2 ±5.3 

4.3 

1.0 x 10 -4 

2.6 x 10 10 

5.6 x 10 +5 

30 

2164, 95% RH 

12.9 ± 1.1 

6 

.06 

2.3 x 10 1 

3.4 xl0 +1 

65 

2164, <2% RH 

22.1 ±3.9 

39 

.01 

3.0 x 10 1 

4.4xl0 +3 

48 

8070 SB, 95% RH 

19.8 ±2.6 

60 

1.6 

8.2 x 10‘ 2 

1.3 xl0 +1 

25 

8070 SB, <2% RH 

25.0 ±3.9 

3,079 

93 

5.5 x 10' 3 

1.5 x 10 +0 

25 

8330, 95% RH 

17.1 ± 1.3 

5 

0.7 

5.6 x 10 1 

1.0 xl0 +1 

25 

8330, <3% RH 

24.5 ±3.9 

19 

0.4 

8.0 x 10' 1 

3.1 xl0 +2 

30 

8330, <1% RH 

30.0 ±3.6 

2,855 

266 

2.1 x 10' 2 

2.3 xlO 10 

30 

Sb-doped. D ' stlllcd 
1 Water 

12.2 ±0.4 

12 

1 

1.1 x 10 1 

1.3 x 10 +0 

115 

Sb-doped, 65%RH 

13.1 ±0.4 

92 

9 

1.7 x 10' 2 

2.4 x 10' 1 

115 

Soda-lime silicate, 
Distilled Water 

20.0 ±2.0 

18 

423 

9.0 x 10' 2 

5.7 x 10 +0 

30 

































































































Table 3 Comparison of fitting ranges and methods for the 8330 borofloat glass tested in 95% RH. 


Fit Method 

n 

B 

MPa 2 s 

A 

m/s (MPaVm)" 11 

# Observations 

All Data (high rate included) 

Individual points 

19.2 ± 1.3 

1 

3.7 x 10° 

30 

Median values 

21.6 ± 3.0 

0.4 

2.1 x 10 1 

6 

Average values 

19.3 ±2.0 

1 

3.9 x 10° 

6 

< 1000 MPa/s (avoid inert region) 

Individual points 

17.1 ± 1.3 

5 

5.6 x 10" 1 

25 

Median values 

19.9 ±3.6 

1 

4.5 x 10° 

5 

Average values 

17.2 ± 1.8 

5 

5.9 x 10 4 

5 


Table 4 Comparison of fitting ranges and methods for the 8330 borofloat glass tested in 1% RH. 


Fit Method 

n 

B 

MPa 2 s 

A 

m/s (MPaA/m)" n 

# Observations 

All Data (high rate included) 

Individual points 

36.8 ±4.4 

608 

5.5 x 10 4 

35 

Median values 

38.2 ± 11.2 

541 

9.0 x 10 4 

7 

Average values 

36.8 ± 8.8 

632 

5.4 x 10 4 

7 

< 1000 MPa/s (avoid inert region) 

Individual points 

30.0 ±3.6 

2855 

2.1 x 10‘ 2 

30 

Median values 

30.4 ± 8.7 

3032 

2.2 x 10‘ 2 

6 

Average values 

30.0 ±7.0 

2984 

2.0 x IQ’ 2 

6 

Table 5 Comparison of propag 

;ation-of-errors equations to Monte Carlo estimates 

Material and 
Humidity 

B 

MPa 2 s 

B-95%o 

MPa 2 s 

A 

m/s (MPaA/m)" n 

A +9S% 

m/s (MPaA/m)" n 

2164, 95% RH 

5.7 

0.06 

0.230 

34.1 

Monte Carlo 

6.1 

0.11 

0.217 

20.2 

2164, <2% RH 

39 

0.009 

0.298 

4,361 

Monte Carlo 

41 

0.011 

0.263 

12,891 


Discussion 


Effects of Humidity 

Table 2 demonstrates that lower test humidity results in higher estimates of n and 
B, regardless of the type of glass tested, implying that controlling or eliminating moisture 
via coatings, etc. will improve component life. The variances are also somewhat larger 
for dry conditions because the shallower slope is more difficult to characterize for the 
same x-axis range. The parameters are also very sensitive to small changes in humidity 
at low humidity: the value of B changes by a factor of ~ 100 for a change of 3% to 1% 
whereas a change from 95% to 3% only results in a factor of -TO change. 







Effect of Fit Method and Range 


The effects of fit range and method on the estimated parameters can be seen in 
Tables 3 and 4: the fitting methods produce similar results for data set; and the inclusion 
of the high stress rate data (1000 MPa/s) substantially alters the results at low humidity 
by increasing n. The lack of an effect of fit method implies either few outliers or 
sufficient data to mitigate the influence of outliers. The effect of fit range can be 
mitigated by using crack growth data only from lower stress rates (< -200 MPa/s) [9], 
and independently measuring inert strength with ~0%RH. 

Confidence Intervals 

The 95% confidence intervals on B for the sealing glasses in Table 2 differ from 
the estimates by 1 to 3 orders of magnitude, even for data sets with 60 observations. The 
relatively large confidence intervals on some of the data sets imply that the use of 
inherent or natural flaws requires very large data sets. In another study, the use of 115 
tests to characterize the Sb-doped glass in Table 1 resulted in standard deviations on n of 
less than % of those for the sealing glasses, and a 95% confident interval on A and B 
within one order of magnitude of the estimate, as shown in Table 2. Improvements can 
also be made by maximizing the range of rates used, and by perfonning most of tests at 
the highest and lowest rates. Monte Carlo estimates compared well with estimates from 
equations (5) to (17), as shown on Table 5. 

Crack Velocity 

The crack velocity as a function of stress intensity based on the estimated 
parameters in Table 2 is shown in Figures 6 and 7. The 8070 glass exhibited the least 
crack velocity whereas the 0120 and 2164 glasses exhibited the greatest velocities at any 
stress intensity. It should be noted that the fracture toughness of some of the glasses was 
not measured and the results could shift if the fracture toughness is significantly different 
than the assumed value of 0.75 MPaVm. Application of common time-to -failure 
equations [7] indicate the sustainable stress for the 8070 is doubled if the humidity is 
changed from 95% to 2%. As compared to soda-lime silicate float glass, the sealing 
glasses exhibit greater susceptibility to slow crack growth, as shown in Figure 6. 

Conclusions 

SCG parameters measured using constant stress rate testing in high humidity are 
systematically smaller (n and B) than those measured in low humidity. Velocities for dry 
environments are - 1 OOx lower than for wet environments: keeping components dry 
should extend the life! The use of high stress-rate data increased estimates n at low RH, 
and thus should be used with caution. 8070 glass-ceramic exhibits the lowest crack 
velocities of the sealing glasses tested, and reducing RH from 95% to 2% nearly doubles 
the sustainable stress. The crack velocity is very sensitive to small changes in RH at low 
RH: for the 8330 glass, a -lOOx change in velocity resulted for changes of RH from 1% 
to 3% and 3% to 95%. Monte Carlo simulations and propagation of errors solutions gave 



similar estimates of parameter variance. Future work should include measurement of the 
fracture toughness of all the materials in order to more accurately estimate A, and 
measurement of the parameters with macro-crack test specimens for comparison. 


95% RH 



Figure 6 Crack velocity for 95% relative humidity based on the parameters in Table 2. 


~2% RH 



Figure 7 Crack velocity for 1 to 3% relative humidity based on the parameters in Table 2. 
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